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Abstract. Purely numerical methods do not always provide an accurate way 
to find all the global solutions to nonlinear ODE on infinite intervals. For 
example, finite-difference methods fail to capture the asymptotic behavior of 
solutions, which might be critical for ensuring global existence. We first show, 
by way of a detailed example, how asymptotic information alone provides 
significant insight into the structure of global solutions to a nonlinear ODE. 
Then we propose a method for providing this missing asymptotic data to a 
numerical solver, and show how the combined approach provides more detailed 
results than either method alone. 



1. Introduction 

Finding global solutions to nonlinear ordinary differential equations on an infi- 
nite interval can be rather difficult. Numerical approximations can be particularly 
misleading, especially because they examine only a finite-dimensional portion of the 
infinite-dimensional space in which solutions lie. Additionally, the conditions for 
global existence can be rather delicate, which a numerical solver may have difficulty 
rigorously checking. In situations where there is well-defined asymptotic behavior 
for global solutions, it is possible to exploit the asymptotic information to answer 
questions about global existence and uniqueness of solutions directly. Additionally, 
more detailed information may be provided by using the asymptotic behavior to 
install artificial boundary conditions for use in a numerical solver. The numerical 
solver can then run on a bounded interval with boundary conditions that match 
the numerical approximations to an asymptotic expansion valid on the rest of the 
solution interval. 

For concreteness, we consider the behavior of solutions satisfying the differential 
equation 

(1) = f"{x) ~ fix) + (j3{x), for all xeR. 

In particular, we wish to know how many solutions there are for a given (j). (There 
may be uncountably many solutions, as in the case where </) = const > 0.) This 
problem depends rather strongly on the asymptotic behavior of solutions to ((T|) as 
|x| — > cxD, so it is useful to study instead the pair of initial value problems 

fo = fix) - p{x) + 0(x) for a; > 
^' \(/(0),/'(0))eZ, 
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and 

fo - fix) - fix) + 0(x) for a; < 
^' \(/(0),/'(0))eZ', 

where (p G C°°(R). The sets Z,Z' supply the initial conditions for which solutions 
exist to ([2]) for all a; > and to Q for all x < 0, respectively. Solutions to ^ 
will occur exactly when Z (1 Z' is nonempty. Indeed, the theorem on existence and 
uniqueness for ODE gives a bijection between points in ZnZ' and solutions to llj. 
[3] Since ([2|) and ^ are related by reflection across a; = 0, it is sufficient to study 
@ only. 

Due to the asympotic behavior of solutions to the methods we employ here 
will be most effective in the specific cases where (j) is nonnegative and montonically 
decreasing to zero. (We denote the space of smooth functions that decay to zero as 
Cjf'(R).) The decay condition on (j> allows the differential operators in ([T]) through 
([3]) to be examined with a perturbative approach as x becomes large, and makes 
sense if one is looking for smooth solutions in L^'(K) with bounded derivatives. 

When <j) is strictly negative, it happens that no solutions exist to ([2]) for all 
X > 0. The monotonicity restriction on (j) provides some technical simplifications 
and sharpens the results that we obtain. This leads us to restrict (/> to a class 
of functions that captures this monotonicity restriction but allows some flexibility, 
which we shall call the M-shaped functions. 

It is unlikely that we will be able to solve ([2]) explicitly for arbitrary 0, so 
one might think that numerical approximations might be helpful. However, most 
numerical approximations will not be able to count the number of global solutions 
accurately. For instance, finite-difference methods are typically only useful for 
finding solutions valid on finite intervals of M. This is unfortunately not sufficient, 
since the behavior of solutions to ^ will be shown in Theorem [7] to either tend to 
zero or fail to exist. A typical finite-difference solution that appears to tend to zero 
may in fact not, and as a result fail to be a solution over all x > 0. 

Because of this failure, we need to understand the asymptotic behavior of solu- 
tions to ^ as we take a:; ^ oo. Equivalently, since (j> G Co"(^)' ^^i^ means that we 
should examine solutions with (p small. The driving motivation for this discussion is 
that solutions to = fix) — f ix) + (pix) for (j) small behave much like solutions to 
= fix) — fix). In the latter case, we can completely characterize the solutions 
which exist on intervals like [a;o, oo). 

In Section [2] we review what is known about the much simpler case where (j) 
is a constant. Of course, then ^ is autonomous, and the results are standard. 
In Section [3l we establish the existence of solutions which are asymptotic to zero. 
Some of these solutions are computed explicitly using perturbation methods in Sec- 
tion [U where low order approximations are used to gather qualitative information 
about the initial condition sets Z and Z' . In Sections [5] and [6l these qualitative 
observations are made precise. Section [7] applies these observations about Z and 
Z' to give existence and uniqueness results for ([1]). Finally, in Section [51 we use the 
information gathered about Z and Z' to provide artificial boundary conditions to 
a numerical solver on a bounded interval, which sharpens the results from Section 
[T] We exhibit the numerical results for a typical family of (j), showing bifurcations 
in the global solutions to ([T]). 



ASYMPTOTIC-NUMERICAL APPROACH FOR ODE 



3 




2. Review of behavior of solutions to = fix) — P{x) + P 
It will be helpful to review the behavior of 

\/(0),/'(0) given, 

where P is a constant, since varying (j) can be viewed as a perturbation on the case 
4i{x) = P. In particular, we need to compute some estimates for later use. We shall 
typically take P > 0, as there do not exist solutions for all a; if P < 0. 



(4) 



Lemma 1. Suppose f is a solution to the initial value problem ^ with /(O) > yP 
and /'(O) > 0. Then there does not exist an upper bound on f{x), when x > 0. 
Additionally, if P <0, there does not exist an upper bound on f{x). 

Definition 2. The differential equation ^ comes from a Hamiltonian, namely 
HifJ') = lf-lf-fP+lp'/'. 

Lemma 3. Suppose f is a solution to the equation (|4|) on R. All bounded solutions 
lie in the funnel 



(5) M - {(/, f')\H{f, /') >Oandf< VP}. 

Any solution which includes a point outside the closure of M is unbounded, either 
for X > or X < 0. (Note that M is the teardrop- shaped region in Figure{^) 

Proof. • M is a bounded set. Notice that H{f,0) > H{f,f'), or in other 

words within M, 

< i/3 - i/'2 -fp+ 2p3/2 < 1^3 ^jp^ 2 3/2 

3-' 2-' ■' 3 - 3'' •' 3 
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Elementary calculus reveals that this inequality establishes a lower bound 
on /, namely that 

(6) - \/3P < / < \/P 

On the other hand, 



(7) l/'l < y ^P3/2 + 2 ;3 _ 2/P < ^ ^p3/4 

immediately establishes a bound on /'. 

• Af is a funnel, from which solutions neither enter nor leave. This is immedi- 
ate from the fact that H is the Hamiltonian, and the definition of AI simply 
says that /') > 0. This sufhces since solutions to ^ are tangent to 
level curves of H . 

• If (/(0),/'(0)) i M then / is unbounded. Evidently if /(O) > VP and 
/'(O) > 0, then Lemma [T] applies to give that / is unbounded. For the 
remainder, discern two cases. First, suppose /(O) > \/P and /'(O) < 
0. Evidently, iJ(/(0), /'(O)) = i7(/(0), -/'(O)), so it's just a matter of 
verifying that a solution curve transports our solution to the first quadrant. 
But this is immediately clear from the formula for 



/' - ±\IIp - 2/P - 2i/(/(0), /'(O)), 



which gives /' = ±/'(0) when / = /(O). The other case is when i/(/(0), /'(O)) < 
0. Then we show that there is a point (V^P, g) on the same solution curve, 
and then Lemma [1] applies. So we try to satisfy 

lp3/2 _ 1^2 _p3/2^ 2^3/2 ^ i/(/(0), /'(O)) < 

= -2i?(/(0),/'(0))>0, 

which clearly has a solution in g. Finally, if 5 = 0, then /(O) > a/P, so it 
has already been covered above. 

□ 

Lemma 4. // / is a solution to ^ with /(O) > VP, and /'(O) > then there 
exists a C such that lim^^^c f{x) — 00. 

3. Existence of asymptotic solutions for e C^f (R) 

The first collection of results we obtain will make the assumption that tends 
to zero. From this, a number of useful asymptotic results follow. Working in the 
phase plane will be useful for understanding Of course ^ is not autonomous, 
but by adding an additional variable, it becomes so. 

Definition 5. We think of as a vector field V on M'^, defined by the formula 
(8) V{fJ',x)=[^^^^{x) 

Notice that the first coordinate of an integral curve for this vector field solves ©. 

Definition 6. Define H{J, f,x) = \p - \f'^ - f(t>{x) + Icfy^'^ix). Notice that for 
constant cf) = P, this reduces to a Hamiltonian for ([4]). 
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Figure 2. The Regions /, //, ///, and IV of Theorem [3 

Theorem 7. Suppose f is a solution to the problem ([2]) where (j) S C^(IR). // / 
does not tend to zero as x —> oo, then there exists a z such that liiRx^z .f{x) — oo. 
Stated another way, if f solves ^ for all x > 0, then hma;_^oo fi^) = 0. 

Proof. If / does not tend to zero, this means that there is an i? > such that for 
each xq > 0, there is an x > xq so that \f{x)\ > R. But since (j) tends to zero as 
X ^ oo, for any P > we can find an xi > such that for all x > xi, \4>{x)\ < P. 
Choose such a P so that the set M in Lemma [3] associated to dH is contained 
entirely within the strip —R < f < R. We can do this since the set M is bounded, 
and its radius decreases with decreasing P, as shown in ([6]) and (O. But this means 
that there is an X2 > xi such that |/(a;2)| > R- 
Construct the following regions (See Figure [2|): 

/ = {(/,/',a:)|/>Pand /' < 0}, 

// = {(/, /',x)|/>Pand /'>0}, 

/// = {(/, f',x)\f<R}, 

and 

IV = |(/, /', x)\f > and /' > -R and Q/^ - ^/'^ - /P + ^p3/2 < if / < x/P 

The following statements hold: 

• Region / is an antifunnel. Along f = R and /' = 0, solutions must exit. 
Once a solution exits Region /, it cannot reenter. Also, because / > -y/P, 
f" ~ f^ — (j) > f'^ — P > 0, solutions must exit Region / in finite x. 

• Region // is a funnel. Along f = R and /' = 0, solutions enter. Now 
/" = /2 — > /2 — P > and /' > 0, so solutions will increase at an 
increasing rate and so, they are unbounded. 
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• Solutions remain in Region /// for only finite x, after which they must enter 
Region IV. This occurs since / < — VP < 0, and so /' always increases. 
Note that for /' < 0, solutions will enter Region /// along / = —R, and 
for /' > 0, solutions exit along / — —R. 

• Region IV is a funnel. Solutions enter along / = — i? and along /' = 
(note that |/| > \/P in both cases). Along the curve boundary of Region 
IV, we have that 

= f'{P-4')<0, 

so that solutions enter. 

Now suppose {f{x2),f'{x2),X2) G /. After finite x, say at a; = x'2, the solution 
through that point must exit Region /, never to return. Then, there is an X3 > x'2 
such that |/(a;3)| > R. So this solution has either {f {xs) , f {x^) , X3) € J/ or S ///. 
The former gives the conclusion we want, so consider the latter case. The solution 
will only remain in Region /// for finite a;, after which it enters Region IV , say at 
X = x'^. Then there is an X4 > x'^ such that |/(a:^4)| > R- Now the only possible 
location for {/{xi), f'(x4), X4) to be is within Region //, since it must also remain 
in Region IV. As a result, the solution is unbounded by an easy extension of 
Lemma [TJ As x becomes large, (j) tends to zero, so the solution will be asymptotic 
to an unbounded solution of = /" — But Lemma [T] above assures us that such 
a solution is unbounded from above, and Lemma |4] gives that it has an asymptote. 
Hence, our solution must blow up at a finite x. □ 

This result indicates that solutions to ^ which exist for all a: > are rather 
rare. Those which exist for all a: > must tend to zero, and it seems difficult to 
"pin them down." We now apply topological methods, similar to those employed 
in [2], to "capture" the solutions we seek. 

We begin by extending the usual definition of a flow slightly to the case of a 
manifold with boundary. 

Definition 8. Suppose M is a manifold with boundary. A flow domain J is a subset 
of M X M such that if a: e J then = pri{J n M x {a;}) is an interval containing 
0, and if x is in the interior of AI then is in the interior of J^. {pr^ : M x Af R 
is projection onto the first factor) 

Definition 9. A (smooth) flow is a smooth map $ from a flow domain J to a 
manifold with boundary M , satisfying 

• $(0, a;) = a; for all a; £ M and 

• ^{ti +t2,x) = <I>(ti, <I>(t2, a;)) whenever both sides are well-defined. 

Additionally, we assume that flows are maximal in the sense that they cannot be 
written as a restriction of a map from a larger flow domain which satisfies the above 
axioms. We call the curve ■ Jx ^ M defined by ^x{t) = ^{t,x) the integral 
curve through x for 

Definition 10. Suppose <i> : J ^ A/ is a flow on AI and x G dAI. Then the flow 
at X is said to be inward-going (or simply inward) if Jx is an interval of the form 
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[0, a) or [0, a] for some < a < oo. Likewise, the flow at x is outward-going if Jx is 
of the form (a, 0] or [a, 0] for — oo < a < 0. 

Theorem 11. (Antifunnel theorem) Suppose ^ : J ^ M is a flow on M and that 
{A, B} forms a partition of the boundary of M such that the flow of ^ is inward 
along A and outward along B. If every integral curve of $ intersects B in finite 
time (ie. Jx is bounded for each x), then A is diffeomorphic to B. 

Proof. For each x E A, Jx — [0,tx], where tx is the time which the integral curve 
through X intersects B. (We have that ^{tx, x) is outward-going, since Jx is closed, 
so it is in B.) 

Using this, we can define a map F : A ^ B hy F{x) = (^{tx^x). Claim that 
F takes A smoothly and injectively into B. The smoothness follows from the 
smoothness of $ and that dM is a smooth submanifold. To see the injectivity, 
suppose F{x) = F{y) for some x,y G A, so ^{tx,x) = ^{ty,y). Without loss of 
generality, suppose < tx < ty. Then we have that 

F{x) = F{y) 

<i>i-tx,F{x)) = ^-tx,F{y)) 

^-tx,^tx,x)) = ^-tx,^ty,y)) 

^{tx-tx,x) = ^{ty-tx,y) 

X = ^{ty-tx,y). 

But the flow is inward at x, so it is also inward at ^{ty — tx,y)- This means that 
(ty — tx — e,y) ^ J for every e > 0. But this contradicts the fact that {ty,y) G J 
unless we have ty <tx- As a result, ty = tx, so x = y. 

In just the same way as for F, we construct a map G : B A so that G takes 
B smoothly and injectively into A. Namely, we suppose Jy = [sy,0] for some Sy, 
and put G{y) = ^{sy, y). Notice that by maximality, if there were to be an a; e A 
such that F{x) = y, Sy = ~tx- 

Now we claim that G is the inverse of F. We have that 

{GoF){x) = ^Sf(x),F{x)) 

= HsF{x),^{tx,x)) 

= ^{SF{x)+tx,x) 

= ^{-tx+tx,x)=X, 

where wc employ the remark about Sy above. □ 

Remark 12. We can extend the Antifunnel theorem to a topological space X on 
which a flow i> : J — > X acts in the obvious way. In that case, there is no reasonable 
definition of the boundary of X . However, the notion of inward- and outward-going 
points still makes sense. If we let A be the set of inward- going points and B be the 
set of outward-going points in X, then the conclusion is that A is homeomorphic 
to B. 

Now wc employ the Antifunnel theorem to deduce the existence of a bounded 
solution to = /" — + (i) for x > xq for some xq > 0. 

Theorem 13. Suppose < <p{x) < K for all x > xq for some > and 
Q < K < oo. Then the region R\ given by Ri = {{f,f',x)\H{f,f',x) > 0,x > 
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Oj / 5: \/ 0(2;)} contains a bounded solution to = f"{x) — f'^{x) + which 
exists for all x greater than some nonnegative Xi . 

Proof. Without loss of generality, we may take = 0, because otherwise solu- 
tions must exit the portions of Ri in {(/, f ^x)\x < xq} since the x-component of 
^(/i f'l ^) is equal to 1. 

If 4>{0) > 0, partition the boundary of i?i into two pieces: A — {{f, f, x)\x — 0} 
and B = {(/, /', x)\H{f, f, x) = 0}. The flow of V is evidently inward along A. As 
for B, notice that Vi? is an inward-pointing vector field normal to B. We compute 

/ - H^) 

\(-.f + VW)W{^ 
= hf + V^W{x), 




which has the same sign as (f>'{x) when / < 4>{x) in Ri. Finally, we must deal 
with the case where / — ^ <j){x) G B. But in this case, f — from the equation for 
H, so we see that V{^y4>ix), 0, x) — (0, 0, l)"^, so the flow is inward when 4>'{x) < 
and outward when 4''{x) > 0. This means that the portion of the boundary of i?i 
on which the flow is outward is a disjoint union of annuli. On the other hand, the 
portion of the boundary of i?i on which the flow is inward is the disjoint union of 
a disk (namely i?i fl {a; = 0}) and some annuli. Now if (j){0) < 0, then the set A 
above is empty. However, it is obvious that set the inflow portion of the boundary 
is still homeomorphic to the disjoint union of a disk and some annuli. 

We can apply the Antifunnel theorem to conclude that there is a solution which 
does not intersect either the inflow or outflow portions of the boundary. There 
is a lower bound on the x-coordinate of such a solution, since the x-component 
of V{f,f',x) is equal to 1, and the Region i?i lies within the half-space x > 0. 
Therefore, there must exist a solution which enters i?i, and remains inside the 
interior of i?i for all larger x. That such a solution is bounded follows from the fact 
that each constant x cross section of Ri has a radius bounded by the inequalities 
(O and (111), and the fact that (f>{x) < K < 00. □ 



4. Asymptotic series solution 

Theorem [13] ensures the existence of solutions to — f" — f^ + (j> for x sufficiently 
large. However, it does not give any description of the initial condition set Z which 
leads to such solutions, nor does it give a description of the maximal intervals 
of existence. Fortunately, it is relatively easy to construct an asymptotic series for 
solutions to which will provide a partial answer to this concern. In doing so, we 
essentially follow standard procedure, as outlined in [1], for example. However, our 
case is better than the standard situation, because under relatively mild restrictions 
this series converges to a true solution. 

We begin by supposing that our solution has the form 

00 

(9) / = E/'^' 

fc=0 

where we temporarily assume fk+i ^ fk and /o ^ 0, as x — ^ +00. (This assump- 
tion will be verified in Lemma [Til) Substituting ([9]) into ([2]), we get 



ASYMPTOTIC-NUMERICAL APPROACH FOR ODE 



9 



E 

fc=0 



fk ^ E/ f'mfk-1 



m=0 



= /o - /o + (/r - 2/0/1 



E 

k=2 



k-1 



fk ^ 2/o/a; — ^ /m/fe- 



We solve this equation by setting different orders to zero. Namely, 
= /0-/0 



= /{'- 2/0/1 + ^ 



fe-i 







fk ^ 2/o/a; — ^ fmfk- 



ni—1 



The equation for /o is integrable, and therefore easy to solve. (There are two 
families of solutions for /q. We select the nontrivial one, because the other one 
simply results in f{x) ~ — J^, (l){s)ds dt.) The equations for fk are linear and 
can be solved by a reduction of order. Thus formally, the solutions are 



fo 



(10) 



f^ = T^[K+rit-drirjm. 

Jk - J [t-a) Jt p^dp 



■ds dt 



ds dt, 



for d, K constants. Notice that these constants parametrize the set of initial con- 
ditions Z. 

Lemma 14. Suppose f{x) = fk{x) where the fk are given by (fTO|) . If there 

exists an M > 0, an R > 0, and an a > 5 such that 

(11) |^(a,)| < _iL_/ora/Z|a:-d| >i?>0, 

then f{x) is bounded above by the power series 



(12) 



^ 00 

l/(-)l^7r-^E 



Proof. We proceed by induction, and begin by showing that the fi term is appro- 
priately bounded: 



l/i(a;)l 



< 



< 



{x - df 
1 



{x - df 

1 



{x - df 



K + J {t-df 
K + J (t-df 
K 



M 



ds dt 



M 



ds dt 



(2 -I- a){5 - a){x ~ d)°'-^ 
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Now since \x — d\ > R and a > 5, we have that 

1 



l/i(^)l < 
< 



\x-d\^ 



\K\ + 



M 



with 
(13) 



\x-d\^' 
Ai = \K\ + 



(2 + a)|5-a|i?°-5 



M 



(2 + a)(a-5)ii«-5 



For the induction hypothesis, we assume that < |^J^|a+» with Ai > and 
for alH < A; — 1. We have that 



k-l 



fc-1 



fmfk-m. < ^ I 



An 



^k—m 



\X - \x - 

^ fe-1 



so by the same calculation as for /i, we obtain 



fk< 



J2m=l A,nAk^ 



1 



(fc + 6)(fc-l) \x-d\ 



fc+2 ■ 



Hence we should take 
(14) 

Hence we have that 



Ak = 



(fc + 6)(fc-l) 



\m\<f:\Mx)\<j^pA, 



fe=0 



X — d 



Lemma 15. The power series given by 



OO A 

frol^-^l" 



with Ao,Ai>{) given, and 



Ak 



X^rn=l AmAk-m _ X^to=1 AmAk- 



□ 



(fc + 6)(fc-l) A;2 + 5A;-6 

converges for \x — d\ > R if Ai < 8R. 

Proof. We show that under the conditions given, the series passes the usual ratio 
test. That is, we wish to show 



lim 



Ak+1 



< R. 
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Proceed by induction. Take as the base case, fc = 1: by the formula for Ak, 

A\ A2 A\ Ax 

Then for the induction step, 

Ak+l _ X]m=l ^mAk-m+1 

~A^ ^ Ak{k^ + 7k) 

_ X]m=2 ^mAk^m+l + AiAk 



< 



< 



< 



< 



Akik^ + 


7k) 


Smil^m+l^fc 


+ AiAk 


Ak{k^ + 


7k) 


^J2m=l ^mAk 


-m + AiAk 


Ak{k^ + 


7k) 


RAk{k^ + 5k- 


6) + AiAk 


Ak{k^ + 


7k) 


R{k^ + 5k-6) 


+ Ai 



(fc2 + 7k) 
i?(fc2 + 5fc-2) 



(fc2 + 7k) 



since Ai < 8R. Thus 



< R for all k, so the power series converges. □ 



Lemma [15] provides conditions for the convergence of the bounding series found 
in Lemma [HI Hence we have actually proven the following: 

Theorem 16. Suppose f(x) = X]feLo/fc(-^) where the fk are given by (jlOp . If there 
exists an M > 0, an R> 0, an a > 5 such that ([TT|) holds, and furthermore 

(15) M < 8(a + 2)(a-5)i?""^ 

then the series for f{x) converges for all x such that \x — d\ > R. 

Proof. Combining Lemmas 1141 and 1151 we find that the key condition is that Ai < 
8R, which by substitution into ([T5|) yields 



< l^^l + 7 ^^/^ r^r. — ^ < 8i?- 

' ' {a + 2){a - 5)R°'-^ " 

But in order to have \K\ > 0, this gives 

M 



0<8R 



(a-|-2)(a-5)i?""5' 

which leads immediately to the condition stated. □ 

Example 17. It is important to notice that the M defined above in Lemma [H] can 
depend crucially upon the value of d and the shape of the curve 4>{x). For the case 
of (j){x) = {x^ — c)e~^ a typical plot of M{d) is shown in Figure[3l It should be 
noted that for various values of c, the M{d) function is numerically very similar. 
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Figure 3. A typical AI{d) function 



1=3erief converges^ □=3erie3 innv diverge c=D.1£, a=1D 




Figure 4. Series convergence test, for (j){x) = (x^ — 0.12)e "^^Z^: 
white = series converges, black = series may diverge 



This also means that the condition p3|) defines a somewhat complicated region 
over which parameters d, K and R yield convergent series solutions. An example 
with our given 4'{x) function is shown in Figure 3) Thus it appears that our series 
solution converges if one goes out far enough, and specifies small enough initial 
conditions. 
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Remark 18. The convergence of the series solution is controned by the convergence 
of a weU-behaved power series. It foUows that as the (p function becomes smaller, 
fewer terms in the series are needed to accurately approximate the solution. Indeed, 
each term in the series solution is asymptotically smaller than the one previous. 
Thus, we can gain some qualitative information from the leading two terms of the 
series, which are 



6 



+ 



1 



(a;-d)2 (x-d)3 
Taking a derivative by x gives 
-12 -3 



{x - d)3 {x - dY 



K 



K 



{t~df 



is - dy 



-ds dt 



[s - d)- 



;ds dt 



+ {x-df 



{s ~ dy 



:ds. 



On the other hand, using the standard expansion for (a + 6)3/2, obtains 



63 



3-62 



[X - d)6 
63/2 



K 



(X - d)3 
63/2 



{x^dy 

jx-df 3-62 
2-63/2 (x-dy 
3.18 ^ 



{t-df 



:ds dt 



K 



{x-dy e^^x-dy 



it 



it-d) 
df 



(s - dy 

6 r '^(^) 



t (s-dy 



1/2 



ds dt 



is~dy 



rds dt 



which leads to 



(16) 



/3/2 + _ rf)3 



(s - dy 



-ds. 



Notice that this equation depends only on d, not K . So from this we should expect 
that the initial data for solutions to be confined to a thin region in the plane x = 0. 
This will be confirmed in Theorem [22] 

Additionally, the relation /' = — ■\/2/3/3/2 holds exactly for the bounded solu- 
tions of = /" - /2. Indeed, in that case, the set Z is {(/, /')|3/'2 = 2/3, /' < 0}. 
So (|16p indicates that the presence of will deflect the set Z largely in the /' 
direction. This is exactly what we show in Section [6l 

5. Restriction to 4> nonnegative and monotonically decreasing 

We now examine what stronger results can be obtained by requiring (\>(x) > 
and (p'ix) < for all x > 0. This can be expected to provide stronger results, in 
particular because the region Ri employed in Theorem [13] acquires a simpler inflow 
and outflow structure on the boundary, and in particular, solutions will exist for 
all X > 0. A collection of four results indicate that all bounded solutions to ^ lie 
within a narrow region. 

Lemma 19. Suppose (f){x) > and (j)'{x) < for all x > 0. Then the region given 
by Ri — {(/, /', x)\H{f, /', x) > 0, X > 0, / < 4>{x)} contains a bounded solution 
to @. 

Proof. Following the proof of Theorem I13[ we partition the boundary of Ri into 
two pieces: A = {(/,/', x)|x = 0} and B = {(fJ',x)\H{fJ',x) = 0}, noting that 
the flow of V is inward along A. Reviewing the computation in Theorem [TS] the 
flow is outward along all of B. 
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Now we employ the Antifunnel theorem, noting that while A is simply-connected, 
B is not. Hence they cannot be homeomorphic, and so there must be a solution 
that remains inside Ri (which evidently starts on A). But the first coordinate of 
such an integral curve must obviously be bounded, since the x cross-sections of Ri 
form a decreasing sequence of sets, ordered by inclusion, and the cross-section for 
cc = is a bounded set. □ 



Lemma 20. Suppose 4'{x) > and 4>'{x) < for all x > 0. Then the region 
given by i?2 = {{f,f,x)\H{f,f,x) < 0, i/^ - i/'2 > 0,x > 0,/' < 0} contains a 
bounded solution to 



Proof. Partition the boundary of i?2 into two pieces: 



A = {(/, /', 0)1/' < 0} U {(/, /', x)\H{f, /', x)^OJ< f < 0}, 

and 

B = {(/, /', x)\Hif, /', a:) = 0, / > /' < 0}U{(/, /', x)\^f-^f'^ = 0, /' < 0}. 

By the calculation in Theorem [131 the flow along A is inward-going. Additionally, 
the flow along the first connected component of B is outward-going. Finally, we 
put S{f,f\x) — i/"^ — i/'^ and observe that \7S is an inward pointing normal 
vector field to B. We compute 

ff'Vf r 

= m^)<o, 

so the flow along this component of B is outward-going. As a result, we can apply 
the Antifunnel theorem, noting that A is connected, while B is not. Therefore, 
there exists a solution to ([2]) that remains in i?2- Note that there is a lower bound 
on the x-coordinate of this solution, since the x-component of V{f, f',x) is equal 
to 1, and the Region R2 lies within the half-space a; > 0. So this solution must 
enter R2 through A, and then never intersect B. Additionally, notice that such a 
solution will have /' < and / > 0, so it must be bounded. □ 



Lemma 21. Suppose (t>{x) > and ip' {x) < for all a; > 0. The complement of the 
set A ^ RiU R2 consists of solutions which are unbounded, and blow up in finite x. 



Proof. Let the complement of the set A be called C, namely C — {(/,/', > 
0} — A. Now the calculations in Lemmas [19] and [20| show that C is a funnel, in 
that the flow through the entire boundary of C is inward. If does not tend to 
zero, then the argument in the proof of Theorem [7| completes the proof, as there is 
a tubular neighborhood about {/ = /' = 0} with strictly positive radius in which 
solutions in C cannot remain. So without loss of generality, we assume — > 0. 
Define the Region / by 



/', x)\f> VW) and (/' > or H{f, /', a;) > 0)} 
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Region 1 

















Region III / 


r 






Figure 5. The Regions /, //, and /// of Lemma [2T] 

There are two bounding faces of Region /, along which the flow is inward. The first 
is 5*1 = / — ^ 4>{x) = 0, along which 




The second was computed already in the proof of Theorem [T31 Notice that /" = 
— (j){x) > in Region I, so f{x) is concave-up, so solutions which enter Region / 
are unbounded. Using similar reasoning to that of Theorem [7l such solutions blow 
up in finite x. 

Now suppose we have a point {a,a',xo) G C with a' < 0. We claim that for 
some Xi > Xq, the integral curve through this point will cross the f — plane. To 
see this, construct Region // by 

1 



II={ifJ',x)qf 



2-' 3 



< and /' < n C. 



Note that 



so the flow is inward along Region // except along f — (along which it is 
outward). Also note that Region // excludes a tubular neighborhood of the line 
/ = /' = with strictly positive radius. As a result of this, the integral curve 
though {a,a',xo) proceeds at least as far as to allow / < —^(j){x), at which point, 
a finite amount of distance in x takes it to /' = 0. 
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So at that point, the integral curve has entered Region ///, say at x ~ xi, where 
III = {(/, ^)\H{f, a^i) < and / < and /' > 0}. 

The flow is evidently inward along /' = and the curved portion by previous 
calculations, and outward along / = 0. Again, note that the line / = /' = is 
excluded from Region /// by a tubular neighborhood of strictly positive radius, so 
there is an X2 > x\ where the integral curve exits Region /// through / = 0. 

Now, consider a point (0,c',X2) along this integral curve with c' > 0. In this 
case, the flow moves such a point rightward. On the other hand, the left boundary 
of Region / moves leftward, approaching / = 0. So there must be an 2:3 > X2 
such that the integral curve through (0, c', x-i) enters the Region /. Collecting our 
findings, we see that every point in C has an integral curve which passes to Region 
/, and therefore corresponds to a solution which is unbounded, and blows up for 
some finite x. □ 

Theorem 22. Suppose (f){x) > and (j)'{x) < for all x > 0. The set Z of initial 
conditions to ^ that lead to bounded solutions 

(1) lies within A ~ RiU R2 and is 

(2) nonempty, 

(3) closed, 

(4) unbounded, 

(5) connected, and 

(6) simply connected. 

(7) Additionally, the portion of Z corresponding to solutions that enter the 
interior of R2 is a 1-dimensional submanifold of {{f, f , x)\x = 0}. 

Proof. (1) From Lemma [211 all bounded solutions must lie in A. 

(2) That there exist bounded solutions in A is the content of Lemmas [TO] and 

(3) Now, put Ao = An {(/, /', 0)} and Bo = dA- Aq. Observe that from the 
proofs of the previous theorems, the flow of V along Ag is inward, and the 
flow along Bg is outward. Since the last component of V does not vanish, 
the flow of V causes each point of Bq to lie on an integral curve starting 
on Aq. This establishes a homeomorphism J7 from Bq into a subset of Aq. 
In particular, is an open map. Now every solution passing through Bq 
is of course unbounded, so Z = Aq — ri(i?o) is evidently closed (it is the 
complement of an open set). 

(4) Bq clearly has the topology of Kx [0,00), S0 7ri(i?o) = 0. Hence, Tri{il{Bo)) = 
also, but notice that 0(i?o) contains OAq. Suppose Z were a bounded 
set. Then it is contained in some disk D. But dD is homotopic to a loop 
in Aq — Z, which either lies in int(Ao — Z) (in which case the homotopy 
need not move it) or in SAq. But this means that the loop encloses all of 
Z, and so cannot be contractible in 0(Bo), which contradicts the triviality 
of 7ri(f2(i?o)). Hence Z is unbounded. 

(5) We first show that the portion of Z lying in the region R2 satisfies the 
horizontal line test. First, note that a solution starting in Z D R2 cannot 
exit i?2. For one, it cannot enter since i?i is an antifunnel. Secondly, 
it cannot exit into M.^ — (i?i U R2) since solutions there are all nonglobal. 
Suppose that /i(0) > ^(O) > and /{(O) = /^(O) with (/i(0), /{ (0)) and 
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(/2(0), /^(O)) both in Z n i?2. But then 

^(/((^)-/2(^)) = /r(^)-/2'(^) 

= fU^)-.fl{x)>o, 



with equahty only if /i(0) = /2(0). Hence, £{fiix) - f2{x)) > for x > 0, 
again with equahty only if /i(0) = /2(0). Now all solutions which remain in 
R2 are monotonic decreasing and bounded from below, so they must have 
limits. On the other hand, the only possible limit is (0,lim2;^oo \/ (t^ix)), so 
therefore all bounded solutions in R2 must have a common limit. Therefore, 
we must have that /i(0) — /2(0). Now this means that the portion of Z 
in the region R2 can be realized as the graph of a function from the /' 
coordinate to the / coordinate. Therefore, if Z were not connected, at least 
one component of Z would be a bounded subset, which is a contradiction. 

(6) Finally, if Z were not simply connected, the Jordan curve theorem gives 
that there are two (or more) path components to 51 (-Bo) ^ — Z , which 
contradicts the continuity of Q. 

(7) By the connectedness of Z and the horizontal line test in R2, the function 
from the /' coordinate to the / coordinate whose graph is Z n int R2 must 
be continuous. Additionally, by the connectedness of Z and the uniquenss 
of solutions to ODE, this implies that the rest of Z whose solutions enter 
the interior of R2 is also a 1-manifold. 

□ 

Definition 23. It is convenient to define, in addition to the initial condition set Z, 
other sets Zx„ C {(/, /', x)\x — xq} such that any integral curve passing through a 
point in Zx^ exists for all a; > 0. Similarly, one can define Z'^^^. 

Remark 24. If — > as x — > 00, we conjecture that Z acquires the structure of a 
1-manifold with boundary. The series solution (jlOp is not valid at such a boundary 
of Z, since such a solution must remain in i?i and therefore decays quicker than 
the leading coefficient of pU]) . Indeed, by analogy with the case where 0=0, 
the leading term /o of the series solution would vanish, and the solution is then 
asymptotic to — (f){s)dsdt. 

All solutions in the form of the series solution (fTO]) enter R2, so a result of 
this theorem is that one of the two parameters d or K in the series solution is 
superfluous. Since d parametrizes solutions when = 0, we conventionally take 
K — 0. Using this, (fT6|) indicates that a good approximation (as xq — > 00, locally 
near / = /' = 0) to the set Z^g is the set 

{i/(/,/')-o} = {(/,/')i^/' = ^n- 

Remark 25. If(/)^P>0asa;^ 00, then it is not true that Z is a 1-manifold 
(with boundary). Indeed, Z has the structure of a 1-manifold attached to the 
teardrop-shaped set M from Lemma |31 

6. Geometric properties of the initial condition set Z 

Lemma 26. Suppose (j){x) > 0, 4>'{x) < for all x > and — s- as x — > cx). 
Then the set Z intersects {(/,/', a;)|/' = 0}. 
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Proof. First, observe that Z intersects the boundary of i?i in a; = 0, since we have 
by Lemmas [19] and [20] solutions entirely within Ri and its complement. Using 
the fact that Z is connected and the Jordan curve theorem, Z must intersect the 
boundary of Ri in the plane x = 0. This reasoning also applies for each Z^g with 
2^0 ^ 0, so that we can find points in the intersections Z^g H dRi for each xq > 0. 
Also note that for the backwards flow associated to our equation (ie. the flow of 
—V), solutions which enter Ri must exit through the plane x = 0. Hence there 
exists a sequence of points {Fn} C Z with F„ = (/n,/4,0) such that the integral 
curve through Fn passes through G„ = G Zn H dRi for each integer 

n > 0. 

Discern three cases: 

(1) If any Fn arc in Quadrants I or II, then since Z is connected, it must 
intersect {/' = 0}. 

(2) If any F„ are in Quadrant III, observe that the flow across the surface 
5 = {(/, /', x)\y^ = /' < 0} is right-to-left. Thus the integral curve 
must cross into Quadrant II on its way to Gn- Therefore, the set Z cannot 
intersect the surface S, and so it must intersect {/' = 0}. 

(3) Assume all the F„ lie in Quadrant IV. Observe that {Fn} is a closed subset 
of i?i n {x = 0}, which is compact. Hence some subsequence of {Fn} must 
have a limit, say F. Since Z is closed, F G Z. But in the portion of i?i 
lying in the x = plane and in Quadrant IV, we have that 

ax 

and 

ax 

Hence f'n> g'n- But since — > 0, 5'^ — > 0, so F lies on {/' = 0}. 

□ 



Lemma 27. Under the same hypotheses as Lemma[2B[ Z also intersects the half 
plane {/ = 0, /' > 0}. 

Proof. Using Lemma [26l we form a sequence {Fn} C Z such that the integral 
curve through F„ passes through {/' = 0, / > 0, a; = n} for each integer n. (This 
can be done without loss of generality, because if any integral curves pass through 
{/' = 0, / < 0}, then the proof is complete by connectedness of Z.) Note that this 
sequence is entirely contained within Ri by Lemma [21] 
Discern three cases: 

(1) There exists an Fn in either of Quadrants II or HI. The result follows by 
the connectedness of Z. 

(2) There exists Fn in Quadrant IV. This cannot occur unless the integral curve 
through Fn passes through Quadrant HI since the flow along {/' = 0} points 
inward into the portion of Quadrant IV inside i?i . 

(3) Otherwise, we assume {Fn} is entirely contained within Quadrant I. In this 
case, note that 

±f' = f^^^<0. 
dx 
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Hence the /'-coordinate of the integral curve through each F„ is positive 
on the interior of Quadrant I. Hence 

:f/ = /'>o, 

ax 

so /„ < gn- But Qn —> Q sincc — *■ 0, so any hmit point of {Fn\ will have 
/-coordinate equal to zero. By the compactness of RiC] {x = 0} and the 
closedness of Z , this implies that Z intersects {/ = 0, /' > 0}. 

□ 



Lemma 28. Suppose > for all x > 0, (/) ^ as x ^ oo, and that there 
exists an xq > such that for all x > Xq, <f>' {x) < 0. Then the set Z intersects 
{/ = 0,/'>0}. 

Proof. We follow the pattern of proving the existence of an intersection for an open 
interval in x containing xq, and then constructing an a priori estimate for the 
/'-coordinate of this intersection. 

Apply Lemma Wf\ to xq, we have that Z^q intersects {/ — 0,/' > 0}. Let 
(0, /q, Xq) lie in this intersection. Note that 

/ = /' > 



and 



dx 



^f' = f-4> = -^<0 

dx 



when evaluated there. As a result, the integral curve passing through (0,/g,xo) 
must pass through Quadrant II first, say for x £ {xi, xq)- Then evidently, Z^-^ must 
intersect {/ = 0, /' > 0}. 

Now since > between xi and Xq, and [a;i,Xo] is compact, there is an 

open set in M'^ containing the intersection of each Z^ with {/ = 0, /' > 0} for each 
X £ [xi,xo], such that in this open set ^/' < i^T < 0. As a result, /{ > /g. Hence 
the /'-coordinate of the intersection point of Z^ with {/ = 0, /' > 0} is decreasing 
with increasing x. (Since we have f"^ — 4> > —(j), it is decreasing at a rate no faster 
than (j). This implies that this intersection point has /'-coordinate no larger than 
(j){x)dx + f^ at X — 0.) Now since solutions through Z^^ exist for all a; > by 
definition, this suffices to show that Z intersects {/ — 0, f > 0}. □ 

Remark 29. The line of reasoning used in the third case of each of Lemmas [26l 
and [57] (and also in [55]) fails if we try to continue Z much farther. This is due to 
the nonmonotonicity of df /dx in Quadrants II and III. More delicate control of (f> 
must be exercised to say more. 

Calculation 30. Towards the end of the more delicate results mentioned in Remark 
[221 it is useful to know the maximum speed along integral curves on points in the 
region i?i in the /- and /'-directions. By this we mean to compute for fixed x the 
maximum values of 



(17) 



I /'I for the /-direction 

1/^ — (j){x)\ for the /'-direction 
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Figure 6. The region A of Lemma [3T] 

in Ri. The first is easy to maximize: we simply look for the maximum value of /' 
in Ri, which is a maximum of 

for — (j>{x) < f < \J This occurs at / = — and has the value of 
\J%l'i<^^l'^ . For the second part of (fT7|l . it is easy to see that the maximum is 3(/)(a;). 
In summary, 

1 l/'l < for the /-direction 

< 30(.t) for the /'-direction 

on 

Using this calculation, we can impose a stronger bound on the decay of 
and constrain the set Z further. 

Lemma 31. Suppose (j){x) > 0, (p'ix) < ~Dj^(j)^/^{x) for all x > for some 

< fc < 1 and D > 1. Then the set Z is contained within {/ > — fc-\/ 0(0)} and 
intersects each vertical and horizontal line in {f > 0} exactly once, and intersects 
{/' = 0} only once. 

Proof That Z intersects {/ = 0, /' > 0} and {/' = 0, / > 0} at aU follows from 
Lemmas mi and [551 Now consider the region A C Ri shown in Figure [5] and defined 

by 

A^Rin ({/' > 0, / < ky/^)} U {/' < 0, 2/3 < 3/'^}) . 

The boundary segments strictly to the right of the boundary labelled 1 in Figure 
|6] are evidently inflow, so long as > 0. The boundary labelled as 1 in the figure 
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moves with speed 

which is greater than maximum speed in the /-direction given in (|18p . This imphes 
that the boundary moves faster than any solution inside Ri. Hence it is an inflow 
portion of the boundary. On the other hand, the curved segment of the boundary 
to the left has been shown to be outflow, in Lemma fT9l 

We observe that the boundary marked 2 in Figure [S] moves with speed 

dx 

which is strictly faster than the boundary marked 1 in Figure [6l and the boundary 
marked 3 in Figure [6] moves with speed 

dxy \ r ^ ' j ^(j,^/^{x) ^ ' 

noting that /'(— ^0(a;)) — ±-y/8/30'^/''(x) is the value of the maximum /'-coordinate 
of i?i at a given x value. This last speed is greater than the maximum speed in the 
/'-direction given by ^ since (t>'{x) < -D%/60^/''(a;). (Notice that ^/6 < 
since < fc < 1.) 

Since D > 1, this means that both the boundaries marked 2 and 3 in Figure 
ini overtake any solution constrained to be within As a result, every solution 
within the region A must leave it within finite x. But the only way to leave A 
causes a solution to enter R'^ — U R2), so every solution which contains a point 
in A cannot exist for all a; > by Lemma 1211 Therefore, Z is contained within 
(i?i U R2) - A. 

Now consider the region B which is defined by 

B^R.n ({/' > 0, / < 0} U {/' < 0, 2f < 3/'2}) , 

which is simply the region A, with k taken to be zero. The portion of the boundary 
of B lying in the {/ = 0} plane is inflow. We can therefore apply the reasoning 
of the vertical hue test: Suppose (/i,/i,0), (/2,/2iO) G Z with /i = /2 > and 
f{ > /2- Then we have both (at x — 0) 

^(/(-/^) = /f-/l-o 

and 

^(/l-./2)=/{-/^>0, 

which gives that /i — /| > for some open interval about a; = 0. Then, -^{f[~f2) ^ 
0, which implies that in fact ^(/i — /2) > 0. However, since all bounded solutions 
tend to the common limit of zero, we have that this implies /{ = /2 at x = 0. 
(Note that since each solution starts in ZD {Ri — B), we have that neither solution 
can become negative, since that would involve entering _B C A or leaving Ri U -R2-) 
This implies that there is a unique intersection of Z with each vertical line. The 
same reasoning applies in the case of the horizontal line test, as in Theorem [221 O 
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Lemma 32. Suppose 4>{x) > for all x > 0, and that (t)'{x) < -D^^^(j)^/'^{x) for 
all X > Xq > for some < fc < 1 and D > 1. Additionally, suppose that for all 
X e [0,xo], 

(19j Xo — X < ;= , 

/S P3/4 
V 3 

where P = maXx£[o.xQ] 4>{^)- Then the set Z is contained within {/ > —\/4>{0)} and 
intersects each vertical and horizontal line in {f > 0} exactly once, and intersects 
{/' = 0} only once. 

Proof. The set Z^q is constrained to lie within the set {/' > ~k\J 0(a;o)}, by Lemma 
|311 (replacing xq by zero) . Now using the /-direction part of ((T8)) , the smallest /- 
value attained in Z^. is 



\\-<i?l\x)dx-k^^) 



If X < Xo, we have 



\\-'i?l^{x)dx-k^ci>{x^) > J-py\x-xo)-k^cj){xo) 



by ([TO)) . As a result, Z^ C {f > —y^lf(x)} for each x < xq. This additionally 
means that in the backwards flow, the entire portion of Z^ contained in {/ < 0} is 
moving away from the plane {/' ~ 0}, which completes the proof. □ 

Remark 33. The condition that <j)'{x) < C(j)^^^{x) implies 

< c 

-^(j)-^^\x) < Cx + C 

4>{x) < J 

{C"-xf 

for some C" and C" . Notice that this condition is satisfied when the series solution 
converges by Theorem 1161 

7. Solutions on the entire real line 

We now combine the results for ^ and ([3]) to discuss properties of the solutions 
to ([1]). When (j){x) is monotonically decreasing, we have by Lemma [211 that the 

initial condition set stays within i?i U i?2. In particular, Z C {/' < y^0^/^(O)}. If 
we relax the restriction of monotonicity, we obtain a similar result. 

Lemma 34. If f — f{x) is a bounded solution to the initial value problem ^ with 
<?!) e C°° n L°°(M) then /'(O) < y^UWli'^ . 

Proof. Since / is a solution to ([2]), then it must satisfy 

r^f-<p{x)>f~Mu. 
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Now Lemma [3] shows that ah bounded solutions to g" — — \\4>\\oo he in the closure 
of the set M given by 

M = |(5,.9')l^g' - ^5" - gUWoo + \ufl' > 0,5 < ^/Mc 

Since this set M is bounded, we can hnd the maximum value of /', which is 

f'ma. = V873||0I|L/'. □ 

Lemma 35. Consider solutions to H]) on the real line, with (p G D L°^{R). If 
for some —00 < A < B < 00, 

4>{x)dx> ^Uli sup \^ix)\f/' + i sup 

y -i \. xe{-oo,A] xe[B,oo) J 

then no bounded solutions exist. 

Proof. For a solution /, we have that /" = — (f>{x) > ~<j){x). Integrating both 
sides we have ^ 

fiB) - f'{A) > - f cj,{x)dx. 

J A 

By Lemma [Ml bounded solutions on 

, x>B have f{B) < y|(sup,g(_^^^] and 

• onx> A, they have f'{A) < (sup^gf^ \(t){x)\ f/'^ , 
so a necessary condition for there to be a bounded solution is that 

''mdx<\[l[i sup \^(x)\f/^ + i sup 



3 \ xe{-oc,A] xe[B,oo) 



□ 



Corollary 36. A necessary condition for bounded solutions to ([T]) to exist if 4> £ 
n L°°{R) is (I){x)dx > 0. 

Proof. Suppose bounded solutions exist. By the proof of Lemma [55l if we let 

gn = - / 4>{x)dx, 



and 

hn^J\({ sup \cl>{x)\fl^ + { sup 
then < for each positive integer n. But the continuity of hmits gives 

(j)(x)dx — hm < hm ha — 0. 

^ n — >OQ ' n — >oo 

□ 

Definition 37. A function (j) e n L°° (M) wiU be called M-shaped if there exists 
an > such that for all |x| > xq, (t>{x) > and 

• (/) is monotonic increasing for x < —xq and 

• is monotonic decreasing for x > xq. 

Theorem 38. Suppose (j) is a positive M-shaped function, then solutions exist to 
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Figure 7. The sets Z and Z' in Example [40] 

Proof. Observe that by Lemma [28l we have that the set Z intersects {/ — 0, f > 
0}. Additionally, by Theorem[211 we have that Z also lies in R2, which is unbounded 
in Quadrant IV. Likewise, the set Z' (for ([3])) intersects {/ = 0,/' < 0}, and 
becomes unbounded in Quadrant I, so Z n Z' must be nonempty, and at least one 
point in this intersection is in the half-plane {a; = 0, / > 0}. □ 

Theorem 39. Suppose cj) is a positive M-shaped function which additionally sat- 
isfies the decay constraints of Lemma \3S\ for x > and x < seperately, then a 
unique positive solution exists to ([T|). (Note that for x < Q, the inequalities and 
signs in Lemma l^ must be reversed, mutatis mutandis.) 

Proof. By the Theorem [551 there exist solutions to ([IJ, one of which comes from 
the intersection of Z O Z' in the half- plane {x = 0, / > 0}. The vertical- line test 
in Lemma [35] allows one to conclude that the solution which passes through that 
half-plane must continue directly to the region R2 of Lemma [20] without crossing 
the plane {/ = 0}. Thus this solution is strictly positive. 

On the other hand. Lemma [32l indicates that Z may lie only in Quadrants I, II, 
and IV, while the set Z' must lie in Quadrants I, III, and IV. On the other hand, 
the vertical- and horizontal-line tests ensure a unique intersection of Z and Z' in 
Quadrants I and IV, so the solution is unique. □ 

Example 40. We examine the family (pdx) = ce~^ which is M-shaped when 
c > 0. Notice that when c < 0, then the necessary condition of Corollary [36] is not 
met, so solutions do not exist for all a; G M. When c = 0, then the trivial solution 
/ = is the only solution. For c > 0, we examine (jj'dx) = ~xce~^ 1"^ . Figure [40] 
shows the sets Z and Z' for the case when c = 0.05. In particular, one notes that 
there appears to be a unique point of intersection. 
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We find the xo for which larger x satisfy < — 4\/2(/)^/^(x)/(fc\/3): 



/c\/3 



which occurs if x > ^^c^^^, so we may take xq = f^C""^/^. 

By way of example, if we fix xq = 4/3, then k = V6c^/-^. (We enforce < A: < 1 
by taking c small.) Now we must check to see if (fT9|) holds. In this case, we need 
to see if c can be chosen so that xq — x ~ 4/3 — x is bounded above by 



v/873ci/4 

g-16/36 _ yg^l/4g-16/36 



> 



/873ci/4 



which can be made as large as one likes by taking c sufficiently small. Noting that 
this last line is a constant in x completes the bound. Therefore, there is a unique 
positive solution for = /" — + ce~^ with c £ [0, e) for some e > 0. 

Remark 41. Taken together, the results of Corollary [36] and Theorems [38] and 
[39] for M-shaped (/> provide the following story about solutions to the equation 
= ./" - + on the real Hue: 

• If the portion of (j) where it is allowed to be negative is sufficiently negative, 
then no solutions exist, 

• If is positive, then a solution will exist. There is no particular reason to 
believe that this solution will be strictly positive or unique. 

• If the decay in the monotonic portions of is fast enough, there is exactly 
one solution, which is strictly positive. 



8. Numerical examination 

8.1. Computational framework. Notice that the results of Remark [41] are not 
sharp: nothing is said if (j) has a portion which is negative, but still satisfies the 
necessary condition of Corollary [JS] Further, if (j) is positive, but does not satisfy 
the decay rate conditions, nothing is said about the number of global solutions that 
exist. Answers to these questions can be obtained by combining the asymptotic 
information we have collected about the sets Z and Z' with a numerical solver. In 
particular, we can obtain information about the number of global solutions to ([T|) 
for any M-shaped cj). 

Suppose that cj) is an M-shaped function, and that xq is such that (l){x) is mono- 
tonic decreasing for all x > xq and is monotonic increasing for all x < —xq. (If 
decreases fast enough, we can choose Xo so that the series solution converges 
on the complement of {—xq^xq) for sufficiently small initial conditions.) Then we 
have the sets Z'_ and Zxg of initial conditions to ensure existence of solutions 
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Plot of (^(x)-(x^-c)e ^ for various choices of c 
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Figure 8. The function 4'{x;c) for various c values 

on (— oo, — a;o] and [a;o,c») respectively. Then any solution to the boundary value 
problem 

{0 = f"{x) - p{x) + </>(a;) for - xq < x < xq 
ifi-xo),fi-xo))eZ'_,^ 
if{xo)J'{xo))eZ^^ 

extends to a global solution of ([T]). So all one must do is solve (|20)) numerically. 
An easy way to do this is to numerically extend the sets Z'_^^ and Z^g to Z' and 
Z respectively (ie. extend them to the plane a; = 0) and compute Z' n Z. 

In order to analyze ([1]) numerically, it is necessary to make a choice of (j). Evi- 
dently, the numerical results for that particular choice of cannot be expected to 
apply in general. However, a good choice of (/> will suggest features in the solutions 
that are common to a larger class of cf). We shall use 

(21) 0(x;c) = (x2-c)e-"'/2, 

where c is taken to be a fixed parameter. (See Figure [H]) This choice of (j) has the 
following features which make for interesting behavior in solutions to ([1]): 

• (j>{x;c) > for c < 0. In this case, there are solutions to by Theorem 
l38l On the other hand, the decay rate conditions are not met over all of R 
so the uniqueness result of Theorem [39] does not apply. Inded, the decay 
rate conditions are met only for sufhciently large |a;|, but not for \x\ small. 

• If c > is large enough, it should happen that no solutions to ([T]) exist, since 
the necessary condition of Corollary [Ml is not met. Indeed, the integral of 
(f) vanishes when c = 1. 

8.2. Bifurcations in the global solutions. Once computed, the numerical solu- 
tions can then be tabulated conveniently in a bifurcation diagram. That is, consider 
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the set in R'^ given by (c, /(O), /'(O)) for each solution /. Evidently, by existence 
and uniqueness for ordinary differential equations, each solution can be uniquely 
represented by such a point. The results of such a computation are shown in Fig- 
ure [9] In this diagram, the solutions are color-coded by the number of positive 
eigenvalues of — 2/. [J] (It should be noted that the green curve continues for 
c < —1.2, but was stopped for display reasons.) 
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Typical global solution curves 
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Figure 10. Typical global solutions: green are from the positive 
branch, the blue one is taken from the lower branch with /'(O) — 0, 
and the red ones are from the fork arms 



Considering the bifurcation diagram, it appears to indicate that ^ undergoes 
a saddle node bifurcation at approximately c = 0.7706, and a subcritical pitchfork 
bifurcation at c = 0.0501. The results agree with Theorem l38l in that solutions 
do exist when c < 0. The saddle node bifurcation was anticipated by the general 
shape of (/). For c > 0.7706, global solutions do not exist, which was qualitatively 
predicted by Corollary [36l 

However, there are some stranger features of the bifurcation diagram. Most 
prominently, the bifurcation diagram appears simply to end near c = —0.4652, and 
at each branch of the pitchfork ate — 0.0740. It is important to verify that these are 
not numerical or discretization errors. If these ends are to be thought of as valid 

,2 

bifurcations, very likely, — 2/ acquires a zero eigenvalue there. Plotting the 
smallest magnitude eigenvalue gives some credence to this possibility. (See Figure 

mi) 

As another check, one can measure the size of the existence interval for solutions 
to H]), centered at a: = 0. Looking in the (c, /(0))-plane (taking /'(O) = 0), one 
can find the first x such that the solution exceeds a particular value. This is 
shown in Figure [T^l in which one sees the same general shape as in the bifurcation 
diagram. (The jagged nature of the graph along the actual bifurcation diagram is 
due to aliasing.) However, for c < —0.4652, the lower branch clearly continues into 
solutions that exist for only finite x. So the end bifurcation indicates a failure of 
the solutions to H]) to exist for all x. 



ASYMPTOTIC-NUMERICAL APPROACH FOR ODE 29 



Smallest magnitude of eigenvalue of numerical approx to (A-2f) 

0.35 I : 1 : 1 : 1 
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Figure 11. Smallest- magnitude eigenvalue measured along the 
lower branch with /'(O) = 




Figure 12. Estimate of existence interval length 



9. Conclusions 

In this article, an approach for counting and approximating global solutions to 
a nonlinear, nonautonomous differential equation was described that combines as- 
ymptotic and numerical information. The asymptotic information alone is enough 
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to give necessary and sufficient (but not sharp) conditions for solutions to exist, 
and provides a fairly weak uniqueness condition. More importantly, the asymptotic 
approximation can be used to supply enough information to pose a boundary value 
problem on a bounded interval containing a smaller interval where asymptotic ap- 
proximation is not valid. This boundary value problem is well-suited for numerical 
examination, and the combined approach yields much more detailed results than 
either method alone. 
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